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ON THE NON-COMMUTATIVE FRACTIONAL WISHART PROCESS. 


JUAN CARLOS PARDO, JOSE-LUIS PEREZ, AND VICTOR PEREZ-ABREU 


Abstract. We investigate the process of eigenvalues of a fractional Wishart process defined by A^ = 
B*B, where B is the matrix fractional Brownian motion recently studied in Il20l . Using stochastic calculus 
with respect to the Young integral we show that, with probability one, the eigenvalues do not collide at 
any time. When the matrix process B has entries given by independent fractional Brownian motions with 
Hurst parameter H G (1/2,1), we derive a stochastic differential equation in the Malliavin calculus sense 
for the eigenvalues of the corresponding fractional Wishart process. Einally, a functional limit theorem for 
the empirical measure-valued process of eigenvalues of a fractional Wishart process is obtained. The limit 
is characterized and referred to as the non-commutative fractional Wishart process, which constitutes the 
family of fractional dilations of the free Poisson distribution. 

Key words and phrases: Fractional Wishart matrix process, measure valued process. Young 
integral, fractional calculus. 


1. Introduction. 


In this paper, we make a systematic study of the dynamics and the limiting non-commutative dis¬ 
tribution of the eigenvalue process of a fractional Wishart matrix process. More specifically, let H E 
(0, l),n,p > 1 and B = > 0}, 1 < i < p, 1 < j < n} be a set of p x n independent 

one-dimensional fractional Brownian motions with the same Hurst parameter H. That is, each bij is a 
zero mean Gaussian process with covariance 


E 




= -{t 
2 ^ 


2H ^ ^2H 




\2H 




As in [l20l . we introduce {N{t),t > 0), the matrix fractional Brownian motion process with parameter 
H whose components satisfy Nij{t) = bij{t), for t >0. 

A fractional Wishart process is the nonnegative definite n x n matrix process defined by JY = 
N*N, where N* denotes the transpose of some matrix N. Let (Ai(f), A 2 (f),An(t),f > 0) be the n- 
dimensional stochastic process of eigenvalues of X and consider the empirical spectral process of the 


eigenvalues Xi {t) > A 2 (f) > • • • > An (f) > 0 of = n ^X, i.e.. 


( 1 . 1 ) 


i=i 


(n) 


(i)’ 


t > 0. 


Different aspects of the dynamics and asymptotics of this spectral process have been considered by 
several authors in the case H = 1/2 of the classical Wishart process. In this case, for n > 1 fixed, Bru 
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Q considered the dynamies and non-eolliding phenomena of the eigenvalue proeess, proving that the 
speetral proeess is an n-dimensional diffusion given by the system of non-smooth diffusion equations 

(1.2) Aj(f) = Aj(0)-f 2 J \/Aj(s) ■ dv\s) + J + 

where v'^ are independent Brownian motions for i = 1,..., n, and denotes the Ito stochastic integral. 
Moreover, Bru Q also showed that if A^"^(0) > • • • > An”^(0), then a.s. the eigenvalues do not eollide at 
any time, i.e.. 


(1.3) 


P(Ai(t) > • • • > \n(t),\/t > 0) = 1. 


The main tool for proving (11.21) is Ito’s formula for matrix-valued semimartingales, and for (11.31) . a 
Mo Kean type argument in the olassioal stoohastio oaloulus. 

Still in the olassioal ease H = 1/2, for fixed t > 0, the asymptotio distribution of is given by the 
olassioal pioneering work of Marohenko and Pastur, [[T4l . Namely, reoall that the/rep Poisson distribution 
(or Marchenko-Pastur distribution) c > 0 is the probability measure on M+ defined by 


where 

(1.4) 




Uc{dx), 

(1 — c)6o{dx) -f Uc{dx), 


c > 1, 
c < 1, 


Uc{dx) = 
a = 


(x — a){b — x) 
271X 

(l-^/c)/ b = 


^{(a,b)}{x)dx, 

(1 + v~cr- 


It was shown in [[T4]| that is the asymptotio speetral distribution, when t = 1, of the empirioal speetral 
measure as lim„_j.oo ^ = c > 0. For fixed positive f 7^ 1, the asymptotio speetral distribution of 
the empirioal speetral measure as lim„_^oo f = c > 0, is the family (/ic(f),f > 0) of dilations of 
the free Poisson distribution whieh is given by /ic(f) = o where ht{x) = tx (see for instanee 
Cavanal-Dubilard and Guionnet BU and Perez-Abreu and Tudor ll25ll ). 

The asymptotio behavior of the empirioal speetral measure-valued proeess of (11.21) falls into the frame¬ 
work of the study of limiting measure-valued proeesses of interaoting diffusion partieles governed by 
Ito stoohastio diffusion equations with strong interactions, as an eigenvalue proeess of a matrix diffusion 
having the property that the partieles never eollide. The general aim in this framework is to show that the 
empirioal speetral measure proeess eonverges weakly in the spaee of oontinuous probability measure¬ 
valued proeesses to a deterministie law. This general direotion of study was oonsidered by Cepa and 
Lepingle [151|, Chan [0, and Rogers and Shi If27]l . among others, in the ease of some Gaussian matrix 
diffusions, turning out to beeome limiting non-oommutative proeesses as a free Brownian motion. See 
also [fT3ll . [fTTll . ifT^ . [l26ll . and referenoes therein. 
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In this tendency, for the case H = 1/2 of the Wishart process, it was proved in dH and [|^ that the 
empirical spectral process {,t > 0 ); n > 1 } converges weakly in the space of continuous probability 
measure-valued processes to the family > 0) of dilations of the free Poisson distribution. The 

proof of this result is mainly based on an appropriate Ito’s formula for matrix-valued semimartingales 
and large deviations or classical Ito calculus inequalities estimates. 

The aim of the present paper is to make a systematic study of the spectra of the fractional matrix 
Wishart process with Hurst parameter H G (1/2,1), and understand several properties such as the dy¬ 
namics of its eigenvalue process, its noncollision properties, and the limiting family of the corresponding 
empirical spectral measure-valued processes. Since fractional Brownian motion is not a semimartingale, 
the main tools we use are based on the Skorokhod and Young stochastic calculus. Recently, Nualart and 
Perez-Abreu lf20ll considered the dynamics and noncollision property of the eigenvalues of a symmetric 
matrix fractional Brownian motion, and in [|2^ . there was derived the functional limit of the correspond¬ 
ing empirical spectral measure-valued processes, which is the non-commutative fractional Brownian 
motion considered by Nourdin and Taqqu [fT^ . 

The final goal of the present paper is to find the non-commutative limit process of the empirical 
spectral measure-valued processes > 0);n > 1} of the fractional Wishart matrix process of 

Hurst parameter H G (1/2,1), asn goes to infinity. Specifically, in this paper we introduce the non- 
commutative fractional Wishart process of Hurst parameter H G [1/2,1) as the family > 0) of 

fractional dilations of the free Poisson distribution given by jjc ^H{t) = o [h^) where hf = t^^x. 
That is. 


(1.5) 




Uc{t){dx), 

(1 — c)6o{dx) + h'cit){dx), 


c > 1, 
c < 1, 


with /ic,H(0) = (5o- Then, as our main results, we prove the following functional limit theorem for 
the empirical spectral measure-valued processes > 0);n > 1}. Let Pr(M) be the space of 

probability measures on R endowed with the topology of weak convergence and let C (R+, Pr(R)) be the 
space of continuous functions from R 4 . into Pr(R), endowed with the topology of uniform convergence 
on compact intervals of R+. 


Theorem 1. Let H G (1/2,1) and > ■ ■ ■ > Xl^\t) > 0,t > 0) be the eigen¬ 
value process of the scaled fractional Wishart process of Hurst parameter H. Assume that 

converges weakly to do, and that lim„_j.oo f = c > 0. Then the family of empirical spectral measure¬ 
valued processes > 0);n > 1} converges weakly in C(R+,Pr(R)) to the unique continuous 

probability-measure valued function corresponding to the law of the non-commutative fractional Wishart 
process of Hurst parameter H > 0) described in t li.5D . 


The strategy to prove this theorem is as follows, including some results that are important on their own. 
We first consider the dynamics and noncollision of the eigenvalues of a fractional Wishart process. The 
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goal of Section 2 is to derive a stochastic differential equation for the eigenvalues of a fractional Wishart 
process in the framework of the Skorokhod integral with respect to the multivariate fractional Brownian 
motion. For preliminaries on the stochastic calculus with respect to fractional Brownian motion, we 
refer to [fT^ . IfTTl . and lfT9l . We start with results on the first and second derivatives of the eigenvalues 
of a nonnegative definite matrix X = N*N as functions of the entries of the matrix N. A detailed 
consideration of the derivatives of the eigenvalues of an Hermitian matrix X, but in terms of the elements 
of X, is carried out in Anderson et. al |l2l and Tao ll2^ . Then we consider in Theorem [21 a new Ito’s 
formula for the Skorokhod integral of functions related to the growth of the second derivative of the 
eigenvalues of fractional Wishart processes, here denoted by X. In Section 3 we prove the noncollision 
of the eigenvalues of X at any time. We follow closely the proof in the case of the fractional symmetric 
matrix Brownian motion in ll20ll . using stochastic calculus with respect to Young’s integral as well as 
appropriate estimates for the moments of the repulsion force of the eigenvalue processes and the joint 
distribution of the eigenvalues of the fractional Wishart matrix. 

The functional asymptotics of the empirical spectral measure-valued process is considered in Section 
4. We first apply our Ito’s formula to find appropriate expressions for the integrated processes 

n 

= t>o. 

^ • 1 
2 = 1 

Then we prove tightness and the weak convergence of the family of measures > 0);n > 

1} in the space C'(M+, Pr(M)).Finally, we characterize the family of laws > 0) of a non- 

commutative fractional Wishart process of Hurst parameter H in terms of the initial value problem for 
the corresponding Cauchy-Stieltjes transform Gc,h of fic,H (see Proposition H]). 

2. The stochastic differential equation for the eigenvalues. 

2.1. Matrix Calculus and notation. In this section, we present some results on the eigenvalues of a 
nonnegative definite symmetric matrix that will be needed during the course of this paper. The compu- 
tatios are similar to those explained in ||2l for the Hermitian and in GOl for the symmetric case without 
the nonnegative definite assumption. 

We denote by J\fpn the collection of p x n matrices. For a matrix N G Afpn, we use the coordinates 
Nij, with I < i < p, I < j < n, to denote the element on the ith. row and the jth column of N. For 
simplicity, we write N = {Nij). Let N* denote the transpose of N. In order to work with Wishart 
matrices, we define X ;= N*N, which is clearly symmetric and therefore belongs to Hn, the space of 
symmetric n-dimensional matrices. 

Let be the set of orthogonal matrices U such that Ua > 0 for all i, Uij ^ 0 for all z, j, and all 
minors of U have non-zero determinants. We denote by Mp^ the set of matrices N G Mpn such that there 
is a factorization 


A := N*N = UAU*, 
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where A is a diagonal matrix with entries A* = Ka such that Ai > A 2 > • • • > A„ and U e . We also 
denote by the space of symmetric n-dimensional matrices X such that X = N*N and N G 
The matrices in the set will be called very good matrices, and we identify with an open subset 
of Moreover, the complement of has zero Lebesgue measure. 

Let Sn be the set 

Sn = • • • ) ^n) G M” : Ai > A 2 > • • • > A„|. 

For any X e Sn, let be the diagonal matrix such that A A = A^. We consider the mapping T : —)■ 

]^n(n-i )/2 defined as follows 

'T'(TJ\ f Ul2 Uip Un—ln \ 


It is known that T is bijective and smooth, see Q, ll^ . Next, we introduce the mapping T : Sn x 
T{U^^) —)■ Hn^ by T(A, Z) = T~^{Z)A^T~^{Z)* which turns out to be a smooth bijection. We denote 
by $ the inverse of T, i.e., $(X) := (A, T(f/)), and observe that it can be defined as a function of the 
associated very good matrix of X: in other words, we define a function ^ Sn x T{U^^) such 

that $(X) := $(X). As a consequence of these facts, it is clear that A(X) is a smooth function of 

Next, we suppose that X is a smooth function of a parameter 6* G M. Then, 



On the one hand, if we compute the values of the eigenvalues of X in terms of the entries of the matrix 
N, we observe 

p n n P / 

( 2 . 7 ) = UriNsrNslUu = EE U'piXgrf 

s=l r=l 1=1 s=l Vr=l 

On the other hand, if we take 9 = we deduce 



( 2 . 8 ) 


dXrl 

dX^H 





Putting all the pieces together, we obtain 


dX^ 

dNkh 


EE EE Uri^kr'^{l=h}Uli 

r=l 1=1 r=l 1=1 


‘^Uhi ^ UriNkr- 

r=l 


Now, we are interested in computing the second derivative of the eigenvalues. We start with the first term 
of (12.61) with 9 = Nkh- Thus from (12.8L we deduce 


d'^Xrl _ 

mf, - 
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implying 


n n 



For the second term of (I2.6L we use again (12.81) and observe 



n n 


^3 l=l r=l 


Therefore 



Now, let us apply the above computations to the particular case of the fractional Wishart process 
which is defined below. Let us consider a family of independent fractional Brownian motions with Hurst 
parameter H G (1/2,1), B = {{bij{t),t > 0), 1 < i < p, 1 < j < n}. As in ||20]| . we introduce 
> 0), the matrix fractional Brownian motion process with parameter H whose components 
satisfy Nij{t) = hjit), for f > 0. 

Definition 1. Let {N{t),t > 0) be the matrix fractional Brownian motion with parameter H. We call 
fractional Wishart process of order n with parameter H to the process {X (t),t > 0) satisfying X (f) = 
N*{t)N{t),fort > 0. 

Following the previous discussion, for any i G {1,...we deduce that there exists a function 
<|)j : —)■ M, which is C°° in an open subset G C with having Lebesgue measure 0, and such 

that Xft) = for f > 0. Therefore using the fact that Nkh{t) = hkh{t), we have 


(2.9) 



and 


( 2 . 10 ) 
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On the other hand, we note 



where in the last identity we have used the faet that U is an orthogonal matrix and the identity (12.71) . 


2.2. Stochastic calculus for the fractional Brownian motion. In order to deseribe the evolution of the 


eigenvalues of a matrix fraetional Brownian motion we present a modifieation of Theorem 3.1 of [|20]| . 
whieh is a multidimensional version of the Ito formula for the Skorokhod integral, in the ease of funetions 
that are smooth only on a dense open subset of Euelidean space and satisfy some growth requirements. 
More specifically, the modification is related to the growth of the second derivative. 

We refer to the monograph of Nualart lfT9ll for the definition of the Skorokhod integral. For the 
definition of the space for p > 1 and 1 < i < n, we refer to section 2 of If20]| . 

Theorem 2. Suppose is an n-dimensional fractional Brownian motion with Hurst parameter H > 
1/2. Consider a function F : M” ^ M such that: 

(i) There exists an open set (7 C M" such that has zero Lebesgue measure and F is twice 
continuously differentiable in G. 

(ii) |F(a;)| + |^| < (7(1 + \x\^), for some constants (7 > 0 and M > 0 and for all x G (7 and 
i = 1,... ,n. 

(iii) For each i G {1,..., n} and for each s > 0 and q G [1, 2), 



for some constant (7 > 0. 

Then, for each i = and t G [0,T], the process {^(Fj^)l[o,i](s), s G [0,T]} belongs to the 


space 



n 


dF 

dxi 


7 r 

f 


n 



( 2 . 12 ) 


Proof. We observe that the proof of this result employs similar arguments as those used in the proof of 
Theorem 3.1 of |[20]| . with the exception of the argument that verifies that Equation (12.121) is well defined. 














8 


J. C. PARDO, J. L. PEREZ, AND V. PEREZ-ARRET I 


To this end, it is enough to verify that the process Ui{s) = Ip,*] (s) belongs to the space ■ 

Indeed, using conditions (ii) and (iii) we have 


E 


\ui{s)\^^^ ds 


Uo 






Uo 


< oo. 


and 


E 


r-T rT 


\D^^'>Ui{s)\^/^ drds 


- 


d^F, 

1/H 

= E 

/ ^ 

Jo 


ds 


Crr.2 

< —T^. 
- 2 


uo JO 

On the other hand, taking g = 1 in condition (iii), we also have for each i G 


E 


\f‘ 

d'^F 

s^^-^ds 

Jo 

dxf 



< oo, for f > 0. 


given the fact that iT > 1/2. As a consequence, all the terms in (12.121) are well defined. 


□ 


2.3. The SDE governing the eigenvalues of a fractional Wishart process. In this section, we are 
interested in studying the dynamics of the eigenvalues of a fractional Wishart process as governed by a 
stochastic differential equation that depends on the Skorokhod integral. 

Recall from Theorem 7.1.2 of [l8]| that the joint density of the eigenvalues of X{t) on Sn, with respect 
to the Lebesgue measure, satisfies 

(2-13) Cnpfl exp (-^)) n 

j=i ^ \ J J 


Theorem 3. Let H G (1/2,1) and {N{t),t > be a matrix fractional Brownian motion with parameter 
H defined as above. Furthermore, let Nff) be an arbitrary deterministic p x n matrix. For each t > 0, 
let Xi,..., \n be the eigenvalues of the fractional Wishart process of order n, X = N*N. Then, for any 
f > 0 and i = 1,... ,n. 


P 

(2.14) Xi{f) = Aj(0) + EE 

k=\ h=l 


I /Hj ■ 

^{N{s))5hHis) + 2H 
obkh 



Xfs) + Aj(5) \ 

Ai(s) - Aj(s) j 


^2H-l 


ds. 


Proof. Without loss of generality, we may assume that iV(0) = 0. Now let us check that {Xi{t),i = 
I,... ,n} satisfies the conditions of Theorem [2l Using (12.71) and the Cauchy-Schwarz inequality, we 
observe 


'£<s>UN{t)) = J2 

i=l i=l 



< 


E EE^i^rW <»iiA'wii5. 


2 = 1 


A=1 r=l 


where 11 ■ 1 12 denotes the Euclidean norm of the columns of a matrix. In particular, for each 1 < i < n, 
we have 


4.(Af(())| < V^\\N{t)\\l 
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On the other hand, using (12.91) and the Cauehy-Sehwarz inequality, we obtain 

2 


d^i 


dh 


'kh 






. r=l 


implying 


d^i 


db 




kh 


<2||iVW|h 


Therefore, for eaeh 1 < i < n, we have 


|^>,(iV(f))| + 




db 


kh 


<^||]V(()||^ + 2||]V(()|h. 


Now, observe that using (12.1 II) we ean verify 




ablk 


(N(t)) 


<2UUt)+2Y, 


i¥^j 


p n 


sEE + 

k=l h=l V i^j 

Xi{t) + Xj{t) 


Er=i uijmMhijt) + uAt)Uhjh)bkr{t)\^ 

|Ai(f) - Aj(f)| 

Eti Ui,{t)UMbki{t) + Uri{t)Uhj{t)bUt)\‘' 


|Ai(f) - Xj{t)\ 


= 2(p + E 


|Ai(f) - Xj{t)\ 


Henee using the joint density of the eigenvalues given by (I2.13L the previous equation, and Jensen’s 
inequality, we obtain for g G [ 1 , 2 ] 


E 




dbl. 


iN{t)) 


< 2'^E 


P 


E 


Xi{t) + Xj{t) 

|Ai(f) - Aj(f)| 


i¥=j 


|Ai(t) + Aj(t)|‘^ 


4v+a'.™e / 

i=i ^ V / / 


I Xi + Xj 

I 


-dl, 


where 7 denotes Lebesgue measure and is a positive eonstant that only depends on q, p, n. Using 
the ehange of variables A* = Pi, we observe the last integral in the previous inequality is a eonstant 
that only depends on g, p, n, thus 


(2.15) 


E 


^2$,; 


dblh 




11 


< K 

— ^^n,p,qy 


where is a positive constant that only depends on g, p, n. The result now follows using Theorem [2] 
and (12.1 II) . □ 
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Remark 1. Let us consider the case H = 1/2, therefore the process {N(t),t > 0) corresponds to a 
p X n matrix of independent Brownian motions. So we have that Equation l l2.i4D takes the form 

" ” A,(s) + A,(s) 


A.(«) = A.(0) + 2 5^5^ / 

k=l h=l 


/■* 

^{N{s))5hkH{s) + 2H 
OOkh Jo 


p 


^ xvo'i _ 


i¥=j 


Xi{s) Aj(s) 


ds. 




Using the fact that the process — — {N{s)) is adapted to the filtration generated by the process N, then 


db 


kh 


the Skorokhod integral coincides with the ltd integral, implying that the stochastic integral in l l2.i4D 
satisfies 


V n 


EE 

k=\ h=l 


'0 


P n 


(N(s))Sh,(s) 


k=l h=l 


'0 


{N{s))-dbkh{s) = 2 ^fis)-dYfis), 


where “■ ” denotes the ltd integral. Using d2.9l) . E* satisfies 


P ri 

y<’ = EE/ 

k=i h=i 


^J2tibki{s)Uii{s) 


s/Xi{s) 


■ dbkh{s). 


Computing the quadratic variation ofY\ we observe 

p n 


y rt 

1 —1 (,—1 Jo 


‘ (Et, ht{s)Uu(s) f 


Xi{s) 


ds = t, 


k=l h=l 

where in the last equality we used d2.7D and d2.9l) . Hence from Levy’s Characterization Theorem, we 
deduce that E* is a Brownian motion. Therefore Equation d2.i4D takes the form 

' Xfs) + Aj(s) 

where i/* is a Brownian motion for i = 1,... ,n, which is the system ofSDE’s obtained by Bru [|3l|. 


A.(i)-A.(0) + 2/ ^U7j■ d„'(s) + (p + Ea’wZ// 


ds, 


3. No COLLISION OF EIGENVALUES. 

In this section we show that, almost surely, the eigenvalues of the fractional Wishart process do not 
collide. Recall that the fractional Wishart process is defined as X = N*N, where N denotes the ma¬ 
trix fractional Brownian motion with Hurst parameter H > 1/2. Let us assume that 9f(0) is a fixed 
deterministic symmetric matrix. 

The following result follows closely the proof of Theorem 4.1 [|20]| . despite the fact that the fractional 
Wishart process is not Gaussian. For the sake of completeness, we provide its proof. 

Theorem 4. Denote by {(Aj(f))t>o, 1 < i < n} the eigenvalues of the fractional Wishart process 
> 0). Assume that Xiiff) > • • • > A„(0). Then, 

P(Ai(t) > • • • > Xn{t),\/t > 0) = 1. 
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Proof. We first assume that for fixed to > 0, we have Ai(to) > • • • > A„(to)- Since the matrix X is 
symmetric, we can use the Hoffman-Wielandt inequality (see flUl) to deduce the following 

n 2 


y] (a.m - AiW) < y^(Ay(()-x 

*j=i 




s = 


2=1 






k=l 


On the one hand, from the previous identity and applying the Jensen inequality twice, we get for r > 2, 


Xft) - Xfs) 



r/2 




^ ^ \^ki{f')bkj{t') ('S) 

I 

P , 

^ ^ \ bkiif^b}^j{f^ ("S) 


*,i=i 


k=l 
n p 


<rf V ^ X] X] \bki{t)bkj{t) - bki{s)bkj{s)\\ 
i,j=l k=l 

On the other hand, observe 

bkiifybkjif^ (s)^ bkiijj) (ybkjif^ bkji^S^ i^kiif') 

Hence using the independence between bki and bkj and the previous identity, we get 


(3.16) 


E 


\bki{t)bkj{t) - bki{s)bkj{s)\ < Cr\t- 


where Cr is a positive constant that only depends on r. Putting all the pieces together, we deduce that for 
any T > 0 and s,t E [to, T], there exists a constant Cn,p,r,T depending on n, p, r, T, such that 


(3.17) E 


n p 


\rH 


|Ai(f) - Ai(s)|'' <n'’ V ^^^^\\bkiit)bkj{t)-bkiis)bkj{s)\'' < Cn,p,r,T\'^ - s\ 

i,j=l k=l 

Since rH > 1, we deduce that the paths of A* are Holder continuous of order (3 for any (3 < H. 

We now follow the same arguments as those used in Theorem 4.1 of ll20l . We consider the stopping 
time 

r ;= inf jf > to : Xft) = Xj(t) for some i 7 ^ j|- 

Observe that r > a.s., and that on the random interval [to, t) the function log(Ai(t) — Xj{t)), where 

i 7 ^ j, is well defined. Since the paths of Xi are Holder continuous for each (3 < H, we can apply the 
stochastic calculus with respect to Young’s integral and deduce that for any t < t AT, 


(3.18) 


log(Ai(t) - Xj{t)) = log(Ai(to) - Aj(to)) + 


Therefore for 1 — 77 < a < A, we obtain 


'to ^i{^) 




to Aj(s) Xj(s) 

Xi(t) — Ai(to) 


d(Xi(s) - Xj(s)). 


'to 


Xi(to) — Xj(to) ’ 
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where 




r(l-Q!)[ VAi(s)-Aj(s) Ai(to) - Aj(to) 

{Us) - X,{s))-^ - {X.{y) - \,{y))-^ 


H“ (y. 


'^0 


(s - |/)"+i 

^ 1 fXi{s)-Xi{t) fUi{s)-Xi{y) 

= F( 5 ) [ (i-,)i-o +(!-“)/ — 


dy 


(y-s) 


2-0 


dy . 


and 

Jj(s) := A-” 

We elaim that 

(3.19) P \Iij{s)\\Jj{s)\ds < oo, for all t > to,i 7^ 

Before we prove the above identity, we first observe that Holder’s inequality with exponents i,q > 1 
such that 1/i + 1/q = 1, imply 


E 




< E 




III 


E 


|J,(5)r 


i/g 


From (13.171) . we deduce that for any jd G (1 — a, i7), there exists a r.v. G with moments of all orders 
such that 


(3.20) 


\Xi{u) - Xi{s)\ < kn,p,TG\s - u\^, for all iG{l,...,n}, 


for all s, M G [to, t] with t < T, and kn,p,T is a positive constant that only depends on n,p and T. The 
above leads to the estimate 


(3.21) 


E 






^ kn,p,r,qP 

for all g > 1 and for some constant k„,p,r,g > 0. In order to estimate E[| 1*^ (s) |^], we consider the integral 
part in the definition of lij and we denote it by 

(Ai(s) - Aj(s))-^ - {Xi{y) - Xj{y))-^ 


UAs) := 


-dy. 


'to 


(s - i/)“+i 

Thus using the same estimates as in the proof of Theorem 4.1 of [|^ . we obtain that there exists a 
constant if > 0 such that 


||if,,(.)||,<if 2 “||G“||,, 

(3.22) 


|||A,(i/)-A,(|/)r-^IU|||A,(s)-A,(.)| 


-ii 


IP3 


'to 


+ lll^t(2/) - ^i(2/)l ip3lll^i(s) - Aj(s)|^ ip2)(s-2/)“^ " ^dy, 
where 6=1 — a, ) = — + — + —, with pj > 1 for i = 1, 2, 3. We choose a, pi, P 2 and ps such that 


I pi P2 P3 ’ 

a 

a>^, 


,213 2 

P 3 <2, p 2 < max < —, - > , 


a a 
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which is possible by taking I and pi close to 1 and since a < \ < (3. In order to prove identity (13.191) . 
we need to estimate 


E 


|Ai(s) - Aj(s)|- 


Reeall that the joint density of the eigenvalues Ai (s) > 
the Lebesgue measure. Then, 


E 


|Ai(s) - Xj{s)\- 


for r < 2 . 

> An(s) is given by (12.131) and that 7 denotes 


= c 


n,p 


nk 


Sn 


{p-n-l)/2 

3 


\Xi - exp 


A, 




n 

j<k 


By making the change of variable A* = and performing the integration, we observe 


E 


|Ai(s) - Xj{s)\ 


-r <r n o 


-2rH 


for a positive eonstant that depends on p and n. In other words, E[|Ai(s) — Aj(s)|“''] is uniformly bounded 
on the interval [to,T]. Therefore for alH 7 ^ j, |Jij(s)|| Jj (s)|ds < 00 , and so the elaim (13.191) follows. 

Finally, the identity (13.191) implies that P(T < r) = 1, otherwise we would get a eontradiction since 
log(Ai(r) — Aj(r)) = — 00 . Therefore as T goes to 00 , we obtain P(r = cx)) = 1. We obtain the desired 
result by letting to go to zero. □ 


4. Functional limit for the Fractional Wishart Process. 

For a probability measure p and a p-integrable funetion /, we write (/i, f) = J f{x)fi{dx). Hence, 
since the empirical measure is a point measure, we have, for / G C^, that 

1 " 

(4.23) (/iS"\/) = -5^/(AS”n^))- 

n ^^ 

i=l 

Therefore, applying the chain rule to the last equation, we get 

(4.24) /) = /) + - E / f {^ t \ s )) dXr \ s ). 

^ i=i Jo 

In order to eonsider the dynamics of the measure-valued proeess > 0), we prove the following 

result. Reeall that the fraetional Wishart process is defined by 

x(”)(t) = for t > 0, 


where (t) = n ^/^iV(f) and N denotes the matrix fraetional Brownian motion (n x p). 
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Lemma 1. Let > 0) be the empirical measure-valued process of the eigenvalues of the fractional 

Wishart process > 0). Then for f G we have 




n p n 






2=1 k=l h=l 


H 



0 


^p^f\dx)p^f\dy)ds 

X y 


(4.25) 


+ 


2Hp 


n 




2H 


n 



f"{x)xs^^ ^p^f\dx)ds. 


Proof We first observe from Il20]| that we can apply Ito’s formula with respect to the Young integral to 
the eigenvalues of the process and get 


(4.26) 


AS”’(i) = 


p n 


Y>(o) + ^EE 

v k=l h=l ' 


/•* 

p dbkh 


{N^^\s))dhH{s), 


for any f > 0 and f = {1, 2,..., n}. Hence using (14.241) and (14.261) . we obtain 

n p n pf 


1 


(4”’,/) = (4”’./) + 45EEE / m(iv'”’W))|^(]v<”i(i>))<;6MW. 

^ • 1 7 7 wO 


2=1 k=l h=l 


Now we will be interested in replacing the Young integrals by Skorokhod integrals in the above expres¬ 
sion. To this end, we prove that the condition of Proposition 3 of [[Hi is satisfied. We will denote by 
the Malliavin derivative with respect to b^h, for each 1 < k < h < n. 

We will first show that 


t pt 


(4.27) 



'0 JO 

In this direction, we first observe 

ft rt 


(f{<!?,{N^^\s)))^{N^^\s))]\s - r\^^-^drds < oo, P-a.s. 

V (Jbkh J 



0 Jo 


( /($,(iV(")(s)))^(iVW(5)) )|s-r|2^-2drds 
OOkh 


(4.28) 


n{2H - 1) Jo 


r(<h,(iv(")(s))) 

1 /•* 

+ 




db 


kh 


(iV(^)(s))) s^^-^ds 


< 92 $, 


n{2H - 1) Jo 
Now, using (12.9!) and the Cauchy-Schwarz inequality, 

Oukh ^ 


/'(■i>,(Af‘">(s)))7i7G(]V(”)(s)))s“-‘<is. 

^’^kh 


(4.29) 




Y1 ) <^Y1 


. r=l 


r=l 
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implying 


E 




/"($^(iVW(s))) ^(iV(")(s)) s^^-^ds 


db 


kh 


4 

< -iiri 

n 


ft n 

Ll 




r=l 


ir{s) 




ds 


t 


iH 


= -^linioo < OO. 


On the other hand, using inequality (12.151) . we obtain 


E 


^'^kh 


< Wf'Woo / E 

^ t^n.pA ,2H 

< < OO. 




dKk 


(]Vl”>(s)) 


„ 2 H -1 


ds 


thus, we conclude 


^ /'(4.(/V<”>(s)))^(/V<”l(s))s“-Vs 


< CX 3 P-a.s. 


Putting the pieces together, we obtain that (14.271) holds. 

Now, we apply Proposition 3 of Q (see also Proposition 5.2.3 of dUl) in order to express the Young 
integrals that appear in (14.261) in terms of Skorokhod integrals. Therefore 

{4"’./) = {4"’./) + T75YYY / /'(4.(A''”'W))^(A''”'WMMf) 

^ i=l k=l h=l -^0 

n p n rt 


77(277-1) 




EEE 

2=1 k=l h=l 



0 ^0 


77f ( /'($,(ivW(s)))^(iv(-)(s)) )|5-r|2^-2drd5. 

OOkh 


Next, we apply identity (14.281) and deduce 

n p n 




(n) 


1 

^3/2 


Jo 


77 


2=1 fc=l /l=l 
n p n pi 


r 7^ /•! 

iEEE / /"(«’i(^“(*))) 

i=l fc=l /i=l 




^ ^0 ^^kh 


2=1 k=l h=l 


dbkh 

kh 


{N^'^\s))] s^^-^ds 


Since is an orthogonal matrix, we observe from (14.291) and (12.71) that 

P n , o I \ 2 P / ri \ 2 

fc=l h=\ / fc=l Vr-=1 / 










































16 


J. C. PARDO, J. L. PEREZ, AND V. PEREZ-ARRET T 


Therefore from the latter identity and applying (12.111) . we deduee 

n p n pf 




n3/2 

2H 


i=l k=l h=l “^0 

i=i ■'« V m 


Af(>.) + Af(.) 
^ \t\s) - \f\s) 




AH 




= <f‘i"’./) + ;^EEE / /'(■S>i(A'‘”>W))^(A'‘”>W)A6, 

«_-| 7 T r -I V 0 Al/l 




n p n pi 


i=l 


d^i 


-Jo 

i=l k=l h=l 


kh (,'5 j 


H 



0 JR2 


{f\x) - ^^i^^\dx)^i^^\dy)ds 
X y 


2Hp 


n 



f'{x)s^^ ^yi!f'\dx)ds+ 


AH 


n 



f"{x)xs^^ ^y^J^\dx)ds, 


where in the last identity we used the definition of This eompletes the proof. □ 

4.1. Tightness. In this seetion, we prove that the family of measures > 0) : n > 1} is tight. 

For this purpose, we first prove the following auxiliary result. 

Lemma 2. Let (A”(f),..., X^\t); f > 0) the eigenvalues of the fractional Wishart process with 

parameter H G (1/2,1). Then for all s^t G f),T],we have 

4 


(4.30) 


E 


^y]iA!">(«)-A<">wi 


(n), 


2 = 1 


< (0 c\t - 


where C is a positive constant. 

Proof As in the proof of Theorem 4, we ean use the Hoffman-Wielandt inequality (see [IHl) to deduce 

n 2 ’A / V 

E (a.'”'(*) - t\s)) < E E ('>« 

i=l i,j=^ \A:=1 

From similar arguments as those used at the beginning of the proof of Theorem 4, we get that there exists 
a constant C > 0 such that for any r G N 


E 


- 4”’(*)4?(*))] <^\t- 


E 




1 2rHrTi2rH 


and, 
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Therefore using the Cauchy-Sehwarz inequality, Jensen’s inequality twiee, and the previous inequality, 
we obtain that the exists a constant K > 0, that depends only on T, such that 


E 






2 = 1 


This completes the proof. 


< 


<4e 


EIE’w-a 




. 2=1 


2-1 


< E 


OJ=l \fc=l 

4" 


i,j=l Vfc=l 




□ 


Theorem 5. Assume that p := p{n) is such that p/n c, as n ^ oo. Then the family of measures 
> 0 ); n > 1 } is tight. 


Proof It is easily seen using (14.231) that for every D<ti<t 2 <T,n>l and f E C^, 


(4.31) 


(A-l”'./) - < ^E |/(E’fe)) - /(A<">(«,)) 


.(’a) 


(n-) I 




2=1 


On the other hand by (I3.20L we know that for each n > 1, the functions are Holder continuous of 
order (3 < H. Therefore, since /' is bounded and applying the Mean Value Theorem, we deduce 

/(AS”'(r)) - /(a1">(s))| < ll/'IU|AS”'(r) - \P(s) . 

Hence using the above estimate, identity (14.301) . and Jensen’s inequality, we obtain 


E 



< E 

(-E /(A!”’te))-/(A‘”>(«.)) 

)1 

- 

- 

\ 77/ 1 

\ 2=1 

/ _ 


< II/ 1 LE 




2=1 


(4.32) 




where in the last inequality we used our assumption and that Cf^T is a constant that depends on f and T. 

Therefore, by the well known criterion that appears in Q (see Prop. 2.4), we have that the sequence of 
continuous real processes /), f > 0 ); n > 1 } is tight and consequently the sequence of processes 

{(/i|"'\ f > 0); n > 1} is tight on (3'(M+, Pr(M)). □ 
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4.2. Weak convergence of the empirical measure of eigenvalues. In the previous seetion, we proved 
that the family of measures > 0); n > 1} is tight on C(M+, Pr(M)), whieh allows us to prove 

now our main result. Before that, we first characterize the family of laws (yUc,r/(f), t > 0) of fractional 
dilations of a free Poisson distribution in terms of the initial value problem of the corresponding Cauchy- 
Stieltjes transforms. We note that the case H = 1/2 was proved in Corollary 3.1 of flU. 


Proposition 1. The family > 0) is characterized by the property that its Cauchy-Stieltjes 

transform is the unique solution to the initial value problem 
(4.33) 


dG 


c,H 


dt 


{t,z) = 2H 




dG 


c,H 


dz 


{t,z) 


2H-1 


Gc,h{0, z) — 


t > 0, 
z e C+, 


satisfying Gt{z) E for z E C"*" and 


lim ri\Gt{ip)\ < oo, for each f > 0. 

rj^oo 

Proof Recall from Section 1 that the family of fractional dilations of a free Poisson distribution is such 
that for eachf > D, pc,H{t) = whereis the free Poisson distribution and= f^^x. 

Therefore, the Cauchy-Stieltjes transform Gc,h{1, z) of the distribution Pc,h{1) satisfies 


z) — 


bc,H{t){dx) 


X — z 


pbPo{h/^) \dx) 


pfP{dx) 


pl’P{dx) 


X — z 


rp-l-2H _ y +2H / rj. _ — ’ 

JU L /v L ./TO db /4/D 


for all 2 ; with lm(z)f^ 0. The above equality implies that 


G,,H{t,z) = —Gff{zt-^^) = 

where G/’P is the Cauchy-Stieltjes transform of the free Poisson distribution. 

If we assume that /ic,iy(0) = Pc, 1 / 2 ( 0 ), then by Corollary 3.1 in [|4|| (see also Proposition 2.1 of If25]| ). 
we know that Gcy /2 is the unique solution to the initial value problem 


dG 


c,l/2 


dt 


(t, z) — G‘/. i/2(^i z) + 1 — C + 2zGcy/2(t, z) 


dG, 


c,l/2 


dz 


(t,z) 


Gc, 1 / 2 ( 0 , z) — —, 
z 


t > 0 , 

e C+. 


Therefore 


dG 


c,H 


dt 


(t, z) = 


dG, 


c,l /2 / ,2H 


dt 


(P^,z) = 2Hi 


2H-1 ^^c,l/2 

dt 


(t,z] 


= 2H 

= 2H 


^c,l/ 2 (t , z) + [1 — C + 2zGcy/2{t ,z)] 


{t,z)={P^,z) 

2H 9Gcy/2 , 2H 


^c,h(^^ ^) + 1 — c + 2zGc^h(1, z) 


dGc,H 

dz 


dz 


(t,z) 


(t^^z) 


2H-1 


2H-1 
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On the other hand, note that 


Gc,h(0, z) 


fic,H{0){dx) 

X — z 


Finally, the uniqueness of (14.331) follows from Corollary 3.1 of [|4l. 


□ 


Now, we are ready to prove our main result, i.e., that the weak limit of the sequenee of the measure¬ 
valued proeesses > 0); n > 1} eonverges to the family of fraetional dilations of the free Poisson 

distribution > 0}. 


Proof of Theorem\T\ From Theorem [51 we know that the family > 0);n > 1} is relatively 

eompaet. Henee, for our purposes, we take a subsequenee f > 0); £ > 1} and we assume that it 

eonverges weakly to (/i*, f > 0), with = 

Now we will eonsider the Cauehy-Stieltjes transform > 0) of > 0) given by 


G7^{z) 



X — z 


For fixed f > 0, let us eonsider the eigenvalues of the matrix XGi){f). We observe that 

and XGi){fXj have the same distribution, thus and also have 

the same distribution. This implies that 


GT{z) 



^lGd(^dx) _ 1 r _ 1 f^nu.-2H 

t2H t2H^lG 


Now by the result of Marehenko and Pastur |[T4ll we know that converges weakly almost surely to 
the free Poisson distribution whieh has Cauehy-Stieltjes transform given by 


GfP{z) 


-{z -f 1 - c) -f y/(2r -f 1 - c)2 - 4c 

Yz 


Henee we obtain that the Cauehy-Stieltjes transform Gt of Ht is given by 


Gt{z) = lim G'^fz) 
1^00 




-2H. 


1 —(t Z -f 1 — c) -I- Y{t~'^^Z + 1 — c)^ — 4c 


— {z -I- (1 — c)f^^) + Y{z + (1 — c)f^^y — 
2zf^^ 


Using the above identity and making some straightforward eomputations, it is easy to verify that {Gt, t > 
0) satisfies (14.331) . and therefore the family {fit, t > 0) eorresponds to the family of fraetional dilations 
of a free Poisson distribution. 

Therefore, we eonelude that all limits of subsequenees of {fiG\t > 0) eoineide with the family 
{pt,t > 0), with its Cauehy-Stieltjes transform given as the solution to (14.331) . and thus the sequenee 
4>o : n > 1} eonverges weakly to {fit, t > 0). □ 
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Remark 2. A more general version of Theorem 1 can be obtained in the case when /Xq converges weakly 
to a measure fio which is not necessarily 5o- 

This is done with a more detailed analysis of equation (14.251) when applied to the deterministic se¬ 
quence of functions 


fj{x) = 


X — Zi 


Zi e (Q X Q) n C+, 


and using a continuity argument. It can be proven that the Cauchy-Stieltjes transform {Gt, t > 0) of 
> 0) satisfies the integral equation 

1 1 


Gt{z) — 


Tojdx) ^ ^ 


(4.34) 


X — z 


2Hc 



{y — zfi {x — zfi J X — y 


X + y 2H-1 


Ps{dx)ps{dy)ds 



..2H-1 


{x — zf 


Ps{dx). 


After some computations it is not difficult to see that the Cauchy-Stieltjes transform {Gt, ^ > 0) o/ 
{pt, i > 0) is the unique solution to the initial value problem given by 


dG 

—^ ^ {l, z) = Gc,172(^5 ^) + 1 — c + 2 zGc,1/2(1, z) 

G,v.(0,.)= 

.L X - z 


dG 


c,l/2 


dz 


{t,z). 


t > 0 , 

2 e C+. 
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